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ABSTRACT 

Large dynamic range numerical simulations of atomic cooling driven collapse of gas 
in pre-galactic DM haloes with T y - lr ~ 10000 K show that the gas loses 90 % and 
more of its angular momentum before rotational support sets in. In a fraction of these 
haloes where the metallicity is low and UV radiation suppresses H 2 cooling, conditions 
are thus very favourable for the rapid build-up of massive black holes. Depending 
on the progression of metal enrichment, the continued suppression of H2 cooling by 
external and internal UV radiation and the ability to trap the entropy produced by 
the release of gravitational energy, the gas at the centre of the halo is expected to form 
a supermassive star, a stellar-mass black hole accreting at super-Eddington accretion 
rates or a compact star-cluster undergoing collisional run-away of massive stars at 
its centre. In all three cases a massive black hole of initially modest mass finds itself 
at the center of a rapid inflow of gas with inflow rates of > 1 Mq yr _1 . The massive 
black hole will thus grow quickly to a mass of 10 5 to 10 6 Mq until further inflow is 
halted either by consumption of gas by star formation or by the increasing energy 
and momentum feedback from the growing massive black hole. Conditions for the 
formation of massive seed black holes in this way are most favourable in haloes with 
T V ir ~ 15000 K and U v ir ~ 20 kmscc -1 with less massive haloes not allowing collapse 
of gas by atomic cooling and more massive haloes being more prone to fragmentation. 
This should imprint a characteristic mass on the mass spectrum of an early population 
of massive black hole seeds in pre-galactic haloes which will later grow into the observed 
population of supermassive black holes in galactic bulges. 

Key words: Cosmology: theory - large-scale structure - black holes physics - meth- 
ods: numerical 



1 INTRODUCTION 



Ifleesl (|l978l ) in his famous "flow chart" has mapped out 
many possible pathways leading to the formation of a 
massive black hole. The discovery of very luminous high 
redshift quasars by the Sloan Digital Sky Survey at 
redshifts greater than 6 is widely seen as an important 
clue in this respect. The presence of super-massive black 
holes (SMBHs) as massive as 3 x 10 9 at these ea rly times 
|Haiman fc Loebll200ll : lFanll200ll ; iFan et alj|2006h requires 
a rapid and efficient build-up of SMBHs, challenging 
models based on Eddington limited growth of stellar mass 
black holes (e.g. IVolonteri et al. I I2OO3I ; IVolonteri fc Reesl 
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200E; lHaim an 2006; lLodato fc Natarai anl 120061 : IVolonteril 



20071 ). The collapse of gas in pre-galactic haloes with virial 



temperatures close to but above the threshold for atomic 
cooling driven collapse is perhaps the most promising 
route to the fast formation of massive black holes in 
the mass range of 10 4 — 10 6 Mq which then can act as 
seeds for the growth o f SM B Hs (e.g. Haehnelt fc Reesl 



19931: lu memura et al 




19931; lEisenstein fc LoebT 
2000 



Oh fc Haimanl 
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1995 
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Early attempts to simulate the collapse of gas in such 
haloes numerically confirmed analytical arguments that the 
gas should not fragme nt efficiently if neither metal nor H2 
cooling are important l|Bromm fc Loebll2003l ). Recently the 
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Table 1. Basic properties of the three simulations (A, B and C) presented in RH08 along with those of two further simulations (D and 
E): boxsize (comoving /i -1 Mpc), starting redshift, collapse redshift, DM particle mass (Mq), spatial resolution (h~ pc), total mass of 
the halo (Mq), the virial radius (h^ 1 pc), circular velocity (km sec -1 ), virial temperature (K), maximum baryon density in the halo 
(cm -3 ), angular momentum parameter A, and temperature at the core of the halo (K). All units are physical units, unless explicitly 
stated otherwise. 



advent of adaptive mesh refinement (AMR) simulations has 
enabled numerical simulations of atomic hydr ogen driven 
collap se with a much improved dynamic range. W ise et al.l 
(|2007l ) and Regan fc Haehnelt (2008, RH08) used the AMR 
code Enzo dBrvan fc Normanll 19951 . Il997l : iNorman fc Brvanl 
ll999l ; l0'Shea et al.ll20o4 T to perform numerical simulations 
of the collapse of metal-free gas at the centre of dark matter 
haloes with virial temperatures of > 10000 K in a cosmo- 
logical context. Both groups found that the gas collapses 
isothermally, becomes self-gravitating and settles into a 
close to isotherm al density profile within the DM halo. 
IWise et alj (|2007| ) pushed the resolution of their simulations 
to study the evolution of the innermost mass shells at the 
very centre of the halo reaching sub- Jupiter mass scales. 
RH08 instead followed the dynamical evolution of a much 
larger fraction of the gas in the halo until it settled into 
rotational support. RH08 found that between 0.1% and 1% 
of the gas at the centre formed a massive self-gravitating 
disc. The discs are surrounded by the ongoing isothermal 
collapse of the outer mass shells which their simulations 
were not able to follow due to the prohibitively short 
dynamical timescales in the dense regions of the disc at the 
centre of the haloes. The simulations contain a wealth of 
information regarding the loss of angular momentum and 
mass in-fall rates during the isothermal collapse in these 
haloes. We here discuss the implications for the expected 
further dynamical evolution of the gas. 

In [J2] we will summarize the main numerical results 
from RH08, discuss the stability of the self-gravitating discs 
which have formed and critically assess the assumption 
that H2 and metal cooling may be suppressed in (some) 
pre-galactic haloes with virial temperatures close to but 
above the threshold for atomic hydrogen cooling driven 
collapse. In 33 we use the dynamical evolution of the gas 
settling into rotational support in our simulations to make 
an educated guess about the further dynamical evolution 
of the gas in the mass shells outside of the discs which are 
still undergoing collapse at the end of our simulations. In 
Sj4]we discuss the implications for the formation of massive 



seed black holes and for the early growth of supermassive 
black holes. In ij5]we discuss our results and conclusions. 



2 AMR SIMULATIONS OF ATOMIC 
HYDROGEN DRIVEN ISOTHERMAL 
COLLAPSE 

2.1 The numerical setup 

The setup of the numerical simulations used here is de- 
scribed in detail in RH08. We will give here a brief sum- 
mary of the salient properties of the simulations. The sim- 
ulations were performed using the adaptive mesh refine- 
ment (AMR) code Enzo originally developed by Greg Bryan 
and Mike Norman at the University of Illinois (Bryan & 
Norman 1995b, Bryan & Norman 1997, Norman & Bryan 
1999, O'Shea et al. 2004). We used Enzo in its AMR mode 
to simulate the collapse of the gas at the centre of DM 
haloes with virial temperatures in the range 13000 — 32000 
K assuming that metal and H2 cooling were suppressed. 
The three simulations of haloes with virial velocities of 
19, 20 and 29 km sec -1 described in RH08 will form the basis 
of our study here. The haloes were selected from simulations 
with box-sizes of 2.0-5.0 h' 1 Mpc. The haloes were resolved 
with ~ 30000 to ~ 120000 DM particles and the simulations 
were run with the maximum refinement level set to 18 (sim- 
ulation A) and 16 (simulations B & C) corresponding to a 
maximum spatial resolution of ~ 0.01 pc in each simulation. 
Further properties of the simulated haloes are listed in table 
[1] In the table the reader will also find the properties of two 
further simulations with smaller and bigger virial tempera- 
tures/velocities which will be discussed in H4.5I 

2.2 Properties of the Isothermal Collapse 

The gas at the centre of the DM haloes cools to a tempera- 
ture of 7000 — 8000 K due to atomic hydrogen cooling and 
settles isothermally into a close to p oc r -2 profile as shown 
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Figure 1. Left-hand Panel: The density as a function of radius for simulations A, B & C. The profile is close to isothermal (n oc r~ 2 ) 
over several decades in radius. Middle Panel: The radial velocity as a function of enclosed gas mass. Notice the sharp drop in the radial 
velocity at the mass of the disc in each simulation. Right-hand Panel: The mass in-fall rate for each simulation calculated from the 
density and radial velocity. 
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Figure 2. Left-hand Panel: The surface mass density of the disc in simulations A-C. Note the exponential density profile within the 
disc. Middle Panel: The rotational velocity of the gas. Right-hand Panel: The Toomre parameter, Q, for the gas in the discs. A value 
of Q > 1 (marked with the long dashed line) indicates that the disc is gravitationally stable. The tick marks in each figure indicate the 
scale radii of the discs. 



in the left panel of Figure[T]for simulations A-C. During the 
collapse, the gas develops supersonic turbulent motions of 
order the virial velocity of the DM haloes. The gas efficiently 
loses angular momentum and settles towards the centre with 
radial velocities, V r , of 5 — 15kmsec about a factor 3-5 
smaller than the virial velocity of the haloes (middle panel 
of Figure [T}. The corresponding mass accretion rates can be 
estimated as 

M = 4ivr 2 p(r)V r (r), (1) 

where p(r) is the density at radius r. The mass accretion 
rates range from a few tenths to about two solar masses per 
year (right panel of Figure [!}. 

The central regions in each halo settle into rotational 
support and form compact, fat, self-gravitating discs with 
a mass of ~ 2 x f 4 M (simulation A) and ~ f x 10 5 M 
(simulations B and C). Note that there is a significant drop 
in the radial velocity of the gas at the mass shell marking 
the outer boundary of the discs. The peaks in the mass ac- 
cretion rate profiles at mass shells outside of the discs are 
due to clumps of high density gas within the gas surround- 
ing the discs. 

In Figure [5] we summarize the basic properties of the 
discs. The discs in simulations A-C have an exponential sur- 
face mass density profile with scale length of 0.075, 0.2 and 



0.27 pc, respectively (left panel). In the middle panel of Fig- 
ure [5] we present the rotation curves of the discs which show 
peak rotation velocities in the range from 25 — 60 km sec - . 
Note that the discs have an elliptical shape and are some- 
what fatter than expected for an ideal isothermal exponen- 
tial disc due to their unrelaxed dynamical state. 



3 THE FURTHER DYNAMICAL EVOLUTION 
OF THE GAS 

3.1 Stability and further dynamical evolution of 
the discs 

9 

We were able to follow the dynamical evolution of the 
discs for several rotation periods and the discs appeared not 
to fragment during this time. The stability of self-gravitating 
discs is governed by the stabilizing forces of pressure on small 
scales and rotation on larges scales. The gravitational sta- 
bility of discs is norma lly characterized by the Toomre pa- 
rameter (|Toomrelll96i ). 



where c s is the sound speed, k = s/2 (V ro t/r) (I + 
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dlny rot /dlnr) 1/2 (|Binnev fc Tremainelll987l ; IOh fc Haiman] 
|2002| ) is the epicyclic frequency, Kot is the rotational veloc- 
ity, S is the surface mass density and r is the radius. For 
values of Q < 1 the disc is expected to be gravitationally 
unstable. In the right panel of Figure[2]we show the Toomre 
parameter as a function of radius for the three discs in sim- 
ulations A-C. Within a few scale radii the discs all have 
Q > 1. The Toomre criterion thus predicts that the discs are 
marginally stable and do not fragment confirming the visual 
impression from the simulations. The further disc evolution 
should be governed by the viscous time scale. 

In the framework of modelling th e viscous evolution of 
a sel f-gravitating disc as an a disc (|Shakura fc Svunvaevl 
1 19731 ) the viscous time scale can be approximated as t v i sc = 
Q _1 (J//_R) _ ' 3 tdyn, where t^ yn is the dynamical time, (H/R) 
is the ratio of disc height to the radius of the disc and /3 — 2 
for a geometrically thin disc (|Lin fc Pringle|[l987h . The discs 
in our simulations are rather fat with (H/R) ~ 0.2 — 0.3 
and the "effective" beta is probably somewhat smaller than 
2. Initially the a param eter should be of ord er 0.1 for the 
marginally stable discs (|Lin fc Pringlel ll990T l . The viscous 
timescale on which the disc evolves should therefore be a few 
tens of dynamical times. Due to the viscous evolution angu- 
lar momentum will be transported outwards. The discs will 
contract, the rotation velocity will increase and gas pressure 
support will decrease making the discs more violently un- 
stable. The a parameter will thus increase rapidly to order 
unity and the discs will stay fat due to the turbulent motions 
induced by gravitational instabilities. Outward angular mo- 
mentum transport will continue, per haps, supported by the 
bars- in-ba rs mechanism proposed by Shlosman et al.l (1 19891 . 
1 19901 ) (see iBegelman et all (|2006l , 120081 ) for a discussion in 
this context). The viscous timescale will then be locked at a 
few times the ever decreasing dynamical time scale. 



3.2 The suppression of metal and Hz cooling and 
the role of star formation 

The most critical assumption in our numerical simulations 
is that cooling by atomic hydrogen and helium is dominant 
which prevents the gas from cooling below ~ 7000 K. The 
temperature of the gas close to the virial temperature of 
the halo is the reason that no efficient fragmentation has 
taken place in the simulations and will presumably also not 
occur during the further collapse of the outer mass shells. 
The temperature of the gas is close to the virial temperature 
of the discs which is important for their stability. Cooling 
by atomic hydrogen and helium is only dominant when the 
H2 and metal abundances are very small. We will now dis- 
cuss the plausibility of this assumption for H2 and metals 
in turn. 

At the densities of the gas in our simulations the H2 
formation timescale is generally shorter than the dynami- 
cal scale. In the absence of a strong ultra-violet (UV) flux 
H2 should form with a universal abundance of xh 2 ~ 10~ 3 
|Oh fc Haimanll2002l ). Efficient dissociation of H2 is thus 
necessary for our assumption, that atomic cooling is domi- 
nant, to be valid. 

H2 molecules can be dissociated by ultra-violet radia- 
tion, either directly by photons with energies greater than 
14.7 eV or as a result of electronic excitation in the Lyman- 
Werner (LW) bands (11.2 eV < E < 13.6 eV). Omukai 



& Nishi (1999), Glover & Brand (2001), Oh & Haiman 
(2002), Bromm & Loeb (2003) and Dijkstra et al. (2008) 
have all discussed in considerable detail H2 formation and 
dissociation in haloes similar to the ones we have simu- 
lated here. UV photons in the LW bands have substantially 
longer mean free paths and should thus be more impor- 
tant for the dissociation of molecular hydrogen. Efficient 
dissociation of molecular hydrogen r equires a UV flux of 
F21, lw > 1000 (|Bromm fc Loeb|l2003l). where F21 is in units 



of lO^ergs-Vm^Hz^sr 



Diikstra et ail (|2008l ) have 
suggested that there may exist a small fraction of haloes 
close to a neighbouring actively star-forming halo where the 
external UV background reaches this critical flux. However, 
they estimate this fraction to be very small, a few times 
10~ 7 . More important is therefore probably internal UV ra- 
diation produced within the haloes. Omukai & Nishi (1999), 
Glover & Brand (2001) and Oh & Haiman (2002) argue that 
as little as one massive star could in principle be sufficient 
to dissociate the entire molecular hydrogen content of our 
haloes. The gas in the discs have rotational velocities of 
30 — 50 kmsec -1 and the gas should thus sit in sufficiently 
deep potential wells to withstand the feedback effects due 
to the resulting supernovae. 

The second assumption we have made in our simula- 
tions is that the halo is virtually metal-free. It is not clear 
at what metallicity metal cooling would become important 
for our simulations. lOmukai et ail l|2008l ) have used one-zone 
models to study this problem and found that a critical metal- 
licity of Z cr > 5x 10~ 6 Zq leads to a sufficiently soft equation 
of state to allow fragmentation via dust cooling at densities 
of n ~ 10 10 cm -3 , somewhat higher than the highest den- 
sities in our simulations. In the absence of dust the critical 
metallicity rises to Z CT > 3 x 10~ 4 Zq above which fragmen- 
tation can occur at substantially lower densities. It is unclear 
to what extent metals produced by stars would actually mix 
with the gas. If the gas in the haloes were efficiently enriched 
with metals fragmentation is obviously expected to occur. 
However, it is then still unclear how efficient fragmentation 
and star formation would be. Even if the gas in the ma- 
jority of haloes with T v i r > 10000A' fragments and forms 
an ordinary star cluster a small fraction of haloes with the 
right conditions to justify our assumption of cooling to be 
dominated by atomic cooling may be all that is needed to 
produce the seeds for the much smaller number of SMBHs. 



3.3 Further dynamical evolution of the outer 
mass shells 

The discs in our simulations sit at the centres of the on- 
going collapse of the gas in the DM haloes. Unfortunately 
the dynamical timescales in dense regions of the discs have 
become prohibitively short and our simulations are thus not 
able to follow the further dynamical evolution of the outer 
mass shells. As discussed in more detail in RH08 the gas in 
the outer parts of the haloes are neither supported by ther- 
mal pressure nor by rotation. The gas is instead marginally 
supported by turbulent pressure and should thus settle fur- 
ther until it either fragments or reaches rotational support. 
We will now investigate at which radius we should expect 
rotational support to set in for the gas in the outer mass 
shells. 

In the left hand panel of Figure [3] we show an estimate 
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Figure 3. Left-hand Panel: The ratio of the rotational velocity, estimated using the inertia tensor (see RH08 for details), to circular 
velocity, Vc- The ratio peaks at the enclosed mass of the disc. A value of > 1.0 indicates rotational support. Bight-hand Panel: The ratio 
of the initial angular momentum, Zi n it, to the final angular momentum, Zfl na l- 



of the level of rotational support as a function of enclosed 
mass at the end of the simulation. The figure shows the ra- 
tio of an estimate of the rotational velocity based on the 
the inertia tensor of the gas (see RH08 for details) and the 
circular velocity, V c . As expected the ratio peaks at the en- 
closed mass of the rotationally supported discs for each of 
the simulations. Note that in simulation B the ratio peaks at 
an somewhat high value of w 2.4 due to the complex dynam- 
ical interactions of two clumps in a multiple clump system. 
In the following we will take a ratio of 1.5 as indicative of 
rotational support. 

In the right hand panel of Figure [3] we show the ratio of 
the initial angular momentum to the final angular momen- 
tum as a function of enclosed mass. We choose as the initial 
value of the angular momentum in each simulation the value 
at the time when the maximum refinement level reached is 
15 in simulation A and 12 in simulations B & C at which 
time the mass shells in question are sufficiently resolved. 

The loss of angular momentum varies between a factor 
of 10 and a factor of 2000 for simulations A to C. Towards the 
outer mass shells the angular momentum loss by the end of 
the numerical simulation becomes progressively smaller due 
to the increasing dynamical time at larger radii. 

We first calculate the radius the outer mass shells would 
collapse to if angular momentum were conserved. As the gas 
in a given mass shell has been losing angular momentum 
rapidly up to the end of the simulation this is unlikely to 
be the case but should give an upper limit for the radius 
at which rotational support is reached. If rotational support 
is reached for Kot(-Rsup) = 1.5V c (.Rsup) and angular mo- 
mentum is conserved the final radius of a mass shell with 
enclosed mass, M, is related to the radius at the end of the 
simulation, i? a im, as, 



-Rsup — 



V-cot (-Rsim) 



^c,gas (^sim ) x 1.5 



(3) 



Note that at the largest radii the gas has not yet become 
self-gravitating at the end of the simulation so V c , S as(r S i m ) = 



GMgas/fsim is calculated for the gas mass only. When ro- 
tational support is reached the gas is always self-gravitating. 
In the left hand panel of Figure U we compare the radius 
as a function of enclosed mass at the end of the simula- 
tion (black curves) with the radius when rotational support 
is reached as given by equation [3] (red curves). Note that 
for the mass shells within the discs our estimate for the 
radius of rotational support can be larger than the actual 
radius at the end of the simulation due to our assumption 
Kot(-Rsup) = 1.5K:(-R S up) which is exceeded in some parts of 
the disc. 

The gas in the outer mass shells is still rapidly losing an- 
gular momentum at the end of the simulations. We have not 
fully understood the angular momentum loss mechanism yet 
but it appears that shocks due to the supersonic turbulent 
motions lead to the dissipation of energy and a redistribu- 
tion of angular mome ntum which is pr esumably efficiently 
transported outwards (|Wise et al.ll2007h . As we know little 
yet about the details of this process we can only speculate 
about the further angular momentum loss of the gas in the 
outer mass shells. In the right hand panel of FigureUwe plot 
the radius the gas would collapse to if the gas in all mass 
shells loses the same fraction of angular momentum as that 
in the outermost mass shell of the disc which has formed 
in the respective simulation. The radius where rotational 
support sets in is then given by 



(3(M) 



Kot(-Rs 



c,gas(R sim ) 



x 1.5 



X R sin 



(4) 



where /3(M) is the extra fraction of angular momentum that 
must be lost by the outer mass shells so that their angular 
momentum loss is the same as that angular momentum lost 
by the gas which is rotationally supported at the disc radius. 
For mass shells within the disc radius we set /3(M) = 1. 

In order to get a feel for how compact the configura- 
tion of the gas in our simulations is compared to observed 
astrophysical objects we also show the typical masses and 
radii of observed compact nuclear star clusters, globular 
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Figure 4. Left-hand Panel: The black curves show the enclosed gas mass as a function of radius at the end of simulation A-C. The red 
curve show the radii to which the gas will fall if the gas becomes rotationally supported with Vrot(flsup) = 1.5 X Vc(iJsup) and angular 
momentum is conserved. Right-hand Panel: In this case the red curves show the radius to which the gas will fall if we assume the gas 
becomes rotationally supported with Vrot(flsup) = 1.5 X Vc(-Rsup) and the gas loses the same fraction of angular momentum as the gas 
which has settled into rotational support in the disc at the end of th e simulations. The blue diamonds give masses and radii of a sample 
of observed compact nuclear star clusters as given bv lMerritj feOOSh. Th e black squares and triangles give masses and radii inferred for 
globular clusters and dwarf galaxies based on data from lBelokurov et al. 



clusters, and dwarf ga laxies (|Belokurov et al . 2007; Merrittl 
l200Sl ; ISeth et alj|2008h . We assume a mass-to-light ra t io of 5 
when converting the luminosities in lBelokurov et all (|2007l ) 
to masses. 

Due to the efficient loss of angular momentum the gas 
in our simulations settles into rotational support in a sub- 
stantially more compact configuration than the stellar com- 
ponent of nearby dwarf galaxies which are hosted by DM 
haloes with virial velocities similar to the ones we have sim- 
ulated here, albeit obviously at z = instead of z — 15. At 
the same mass the radii are typically one to three orders of 
magnitude smaller. Observed globular cluster and nuclear 
star cluster are more compact and are generally hosted by 
galaxies/DM haloes much more massive than the DM haloes 
we are considering here making a comparison less straight 
forward. The inner mass shells (< 10 5 M Q ) of the gas m our 
DM haloes settle into a configuration which is also substan- 
tially (by an order of magnitude in radius) more compact 
than globular clusters of the same mass. The configuration 
of the gas in the outer mass shells is expected to settle into 
rotational support at radii which are similar or smaller than 
those of compact nuclear star clusters depending on how 
much angular momentum loss will occur. Note, however, 
that the observed cl usters may have expanded due to the 
rapid removal of gas (jBastian fc Goodwinll200"r} ). 



4 IMPLICATIONS FOR THE FORMATION OF 
MASSIVE SEED BLACK HOLES 

4.1 Pathways to a massive black hole 

Various mechanisms by which a massive black hole may 
form have been laid out comprehensively by Rees (|Rees! 
1 19781 . Il984h . Since then considerable effort has gone into 



researching the paths in this flow chart but the actual for- 
mation process of SMBH nevertheless still eludes us. Fig- 
ure summarizes the paths to a massive black hole con- 
centrating on the possibilities that pre-galactic DM haloes 
with T vir > 10000 K offer. We thereby use the results of 
our numerical simulations as a starting point. In the now 
well established hier archical paradigm for galaxy formation 
l|White fc Reedll978h DM haloes with T vir > 10000 K build- 
up through the merging of several DM haloes with smaller 
virial temperatures. In these less massive haloes gas can only 
collapse and form stars, if the gas cools due to Hi and/or 
metals. Especially for the first generation of these haloes it is 
thus very uncertain, how efficiently they have formed star s 
(e.g. iGreif et al.1 |200S| ; iNormanl [20081 ; IWhalen et alll2008h . 
As discussed in i]3.2l this also introduces considerable un- 
certainty for the cooling processes and star formation effi- 
ciency in the more massive haloes with T v j r > lOOOOJf. If 
fragmentation and star formation is efficient early on, an 
ordinary star cluster/dwarf galaxy may form. If cooling on 
the other hand is dominated by atomic cooling the gas is 
expected instead to s ettle into a rotationally supported, fat , 
self-gravitating disc (|Mo et al.|[l99Sl ; lOh fc Haimanl I2002T I. 
The further fate of this disc depends crucially on whether 
atomic cooling remains the dominant cooling process dur- 
ing the further evolution of the disc. If this is the case the 
gas will not efficiently fragment and will continue to con- 
tract on a timescale controlled by the rate at which the 
gas loses angular momentum. As the gas will be gravita- 
tionally unstable this is expected to occur on a timescale 
which is longer by a factor of a few than the dynamical 
timescale. The further fate then depends on how efficiently 
the entropy produced by the release of gravitational en- 
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Figure 5. Summary of the possible pathways to massive black holes via a stellar seed black hole, a quasi star or via a nuclear star cluster 
in DM haloes with T vir > 10000 K. 1 



ergy is trapped. One possibility is that th e gas settles into 
an high entropy super-massive star (e.g. iHovle fc Fowlerl 



1 The four pictures of simulated haloes at the top of the figure 
and the disc show simulations d escribed in RH08, the sketch of 
the quasi-star is based on that in Bcgclman ct al. 1 1120081) and the 
picture of the star cluster shows the globular cluster M80 (source: 
http:/ /hubblesite.org/ncwsccntcr/archive/1999/26/imagc/a ). 



1963al 


bl: IChandrasekharl Il964allbl: IZel'Dovic 


i & Novikovl 


19671: 


Zel'Dovich Il97fj|; IShaDiro & Teukolskvl 


19791. Il983h 



which ma y take the form of a q uasi-star as recently sug- 
gested by iBegelman et al. | (|2006f ). If the released gravita- 
tional energy is instead transported away in a disc-jet con- 
figuration then rapid accretion onto a stellar mass seed black 
hole would be the expected outcome. Finally if either metal 
or H2 cooling become efficient in the self-gravitating disc and 
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Figure 6. Left-hand Panel: The black curves show the rotational velocity as a function of enclosed mass at the end of simulation A-C. 
The red curves show the rotational velocity the gas attains if it settles into rotational support with Viot (-Rfinal) = 1-5 x Vc(-Rfinal) an d 
angular momentum is conserved. The values of i?fl na i in Figure [4] are used. Right-hand Panel: In this case the red curves show the 
rotational velocity of the gas attains if we assumes the gas reaches rotational support with V ro t (i?ft n al) = 1.5 X Vc(-Rfinal) an d loses the 
same fraction of angular momentum as the gas which has settled into rotational support in the disc at the end of the simulations. The 
values of i?,fi na i in the right hand panel of Figure [4] are used in this case. 



fragment ation is efficient a very compact star cluster should 
form (e.g. lOmukai et al.ll2008l : lBernadetta fc Volonteril2008n 
. In the following we will discuss these different possibilities 
in more detail. 



4.2 Super-massive stars and quasi-stars 

At the time when the gas settles into rotational sup- 
port in our simulation the surface mass density at the 
center exceeds 30g/cm 2 (Figure 1) and the gas is thus 
optically thick. The dominant energy production process 
is the release of gravitational energy. As argued in H3.ll 
this will occur on a few times the dynamical time scale. 
The corresponding mass inflow rate and luminosity rise 
rapidly as the gas contracts. They can be estimated as 
M ~ Q.lvl/G ~ l^c/SOkmsec-YMoyr- 1 and L grav ~ 
Q.lv^/G ~ 3.6(i; c /30kmsec" 1 ) 5 x 10 38 ergs _1 . If the gas 
in the discs in our simulations is able to contract by another 
factor 100 in radius without efficient fragmentation the cir- 
cular velocity will rise to 300 km sec -1 . At this stage the 
luminosity due to the release of gravitational energy will be- 
come comparable to the Eddington luminosity Z/Edd, where 
the outward pressure due to the radiation balances the grav- 
itational attraction, 



L Edd = 47rGMmgC ^1.3xl0 43 ( 



M 



10 5 M< 







-)ergs \ 



(5) 



where c is the speed of light, run is the mass of a hydro- 
gen atom and or is the Thomson scattering cross-section. 

If the entropy produced due to the release of gravita- 
tional binding energy can be trapped i n a quasi-spherical 
confi guration the disc should puff up (e.g. IBegelman fc Reea 
1 19781 ) when the luminosity due to the release of gravita- 
tional binding energy exceeds the Eddington luminosity and 
a (radiation) pressure supported super-massive star should 



ensue. IBegelman et al.l ((2006, 2008 ! ) have argued that a su- 
permassive star embedded in a spherical accretion flow with 
mass accretion rates in the range ~ 0.1 — lA^yr" 1 will 
develop a core-envelope structure akin to that of a Thorne- 
Zytkow object where a co mpact neutron star is e mbedded 
in a diffuse red giant star dThorne fc Zvtkowlll977f ). 

iBegel man et al.l (|2006h dubbed such an object a "quasi- 
star" and envisage its evolution as follows. Initially a com- 
pact core of about 10 — 20 Mq surrounded by a more diffuse 
radiation-pressure supported envelope forms. In the cen- 
tral regions nuclear reactions may commence but the corre- 
sponding energy release is not expected to exceed the release 
of gravitational binding energy. The high mass in-fall rate 
contin ues to compress and heat the core. Bcgcl man et al.l 
|200cf ) argue that the core can reach temperatures of up to 
5 x 10 s K at which point the gas should cool catastrophically 
due to thermal neutrino cooling leading to a core-collapse 
black hole. At this stage the black hole is surrounded by a 
massive gas envelope of more than one hundred times the 
mass of the black hole and the central black hole accretes 
at a rate which corresponds to the Eddington rate of the 
mass in the envelope. Given that the mass of the envelope 
is much larger than the mass of the black hole this allows 
the black hole to quickly grow t o a m ass of 10 3 - 1O 4 M , 
at which point IBegelman et al.l (|2008l ) argue the quasi-star 
evaporates as the photospheric temperature falls. Whether 
this is actually the correct picture is still rather uncertain 
but definitely merits further investigation. It is, however, 
interesting to note that the postulated in-fall rates in the 
Begelman et al. model are similar to those in our simula- 
tions. The in-fall rates of the gas in the outer mass shells at 
the end of our simulations range from 0.3 — 1.5 solar masses 
per year and as we have discussed at the start of this section 
the inflow rates at the centre of the discs, which have formed 
in our simulations, are expected to become even larger as 
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Enclosed Mass [Msolar] Enclosed Moss [Msolar] 

Figure 7. Top Panel: Left-hand Panel: The black curves show the expected collision times (equation [6} for stellar collisions in a compact 
star cluster with masses and radii as that of gas at the end of simulation A-C. The red curves are the expected collision times assuming 
that the gas settles into rotational support without losing further angular momentum using the values of i?fl na i from the left hand panel 
of Figure|4] Right-hand Panel: The collision times assuming that further angular momentum is lost as described in section 3.3 and using 
values of i?fi na i as in the right hand panel of Figure [4] Bottom Panel: The same plots except in this case the relaxation time (equation 
O is plotted against enclosed mass. 



the self-gravitating discs evolve on a viscous timescale due 
to gravitational instabilities. 



4.3 Super-Eddington accretion of stellar mass 
black hole seeds 

If some star formation occurs in the discs at the centre of our 
haloes a stellar mass black hole with a mass of up to 100 Mq 
should easily form. Alternatively a stellar mass black hole 
which has formed in the early stages of the collapse or in one 
of the precursors of our DM haloes may sink to the centre 
of the DM potential in which our discs form. Independent of 
whether a supermassive star will form from the gas settling 
to the centre of the disc such black holes find themselves 
in an ideal environment for rapid growth. The surrounding 
material is flowing to the centre at highly super-Eddington 



accretion rates compared to the Eddington accretion of a 
stellar mass black hole. The rotational velocities of our discs 
range from 30 — 60 km sec -1 and the material is thus already 
strongly bound. The energy released by a stellar mass black 
hole will not be able to halt the surrounding inflow of mat- 
ter and the black hole should grow quickly at highly super- 
Eddington rates to masses comparable to the disc masses if 
fragmentation is not important. In Figure[6]we plot the rota- 
tional velocities at which the outer mass shells are expected 
to settle into rotational support for our two assumptions re- 
garding angular momentum loss. If the gas in the outer mass 
shells lose angular momentum as efficiently as the gas in the 
discs in our simulations the rotational velocity could be as 
large as several hundred km sec -1 . Thi is quite remarkable 
considering that the virial velocities of the DM haloes are 
in the range 20 — 30 km sec -1 . There is thus potential for 
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rapid further growth and the final black hole mass would 
probably be set by competition between gas consumption 
by star formation and accretion for the outer mass shells or 
by the increasing feedback due to the increasing accretion 
luminosity. 

If the entropy produced due to the release of gravita- 
tional energy during the accretion is trapped in a quasi- 
spherical configuration then the accreting black hole would 
actually be very similar to the later stages of a supermassive 
star discussed in the last section. The released gravitational 
energy may, however, not efficiently couple to the surround- 
ing inflowing matter. It could instead be transported out- 
wards in the form of a jet as is observed in young stars, ac- 
tive galactic nuclei and probably also gamma ray bursts (e.g. 
Mundt et aDll990l : iBegelman et aI1ll984l ; iMeszaros fc Reesl 



200lh . Accretion would then occur in a fat accretion disc at 



super-Eddington accreti on rates possib l y involving the pho - 
ton bubble instability ijGammiel Il998l : burner et al.1 120051 ) 
and/or advection dominated in- / outflows as discu ssed by 
iBlandford fc Begelmanl (|l999h and lBegelmanl (|200ll ). 



4.4 The formation and evolution of a compact 
star cluster 



As briefly discussed in £13.21 the evolution of the gas in dy- 
namical haloes with virial temperatures T v j r > 10000 K de- 
pends crucially on whether metals and/or H2 are present 
in sufficient abundance to facilitate cooling well below the 
virial temperature. If this is the case the gas at least in 
principle can fragment and form stars, even though it is un- 
clear how efficient frag mentation and star fo rmation in this 
situation would be (e.g. lOmukai et aTll'illllS 1. If star forma- 
tion occurs early on in the collapse the outcome is proba- 
bly similar to nearby dwarf galaxies, which have DM haloes 
of similar virial velocity and temperature as those we have 
found in our simulations. Some of the nearby dwarf galax- 
ies should have actually formed at similarly high redshifts 
(2 ~ 15) as we are discussing here. Early fragmentation 
requires the presence of a substantial metallicity and most 
likely also dust. Dwarf galaxies have typical metallicities of a 
hundredth to a tenth of solar. Star formation appears to be 
very inefficient in dwarf galaxies. This is generally assumed 
to be due to feedback from supernovae and or ionizing ra- 
diation which appears to strongly affect star formation in 
haloes with V c < 50 km sec . 

It seems, however, likely that at z ~ 15 (and perhaps 
even at the present day) a fraction of DM haloes with virial 
temperatures T v i r > 10000 K have not yet been significantly 
polluted with metals or that the metals in them have not 
been sufficiently mixed into the cold gas phase for metal 
cooling to be relevant. 

In such haloes stars could probably only form by H2 
cooling once rotational support has set in. As discussed in 
S14.ll this would lead to very compact star clusters which 
may well be substantially more compact (see Figure [4} than 
the globular clusters and the compact nuclear star clusters 
which have been discovered in the centres of many galaxies 
including luminous elliptical galaxies and the bulges of spi- 
ral galaxies thanks to the unrivaled resolution of the Hubble 
Space Telescope (e.g. ICarollo et ai1ll997l . 1 19981 : iBoker et all 



120021 : ICote et al.ll2006t l 



star cluster, after rotational support sets in, then the typi- 
cal velocity dispersions of the star cluster would be substan- 
tially larger than the velocity dispersion of the DM haloes 
which host them and could in principle reach several hun- 
dred kmsec -1 if the outer mass shells lose angular momen- 
tum as efficiently as the gas in the discs in our simulations. 

Runaway collision at the the centre of dense stellar 
systems is another promising pathway to the formation of 
a massive black hole which has been studi ed intensively 
(Begelman & Rees 1978; see iFreitad (|2008h for a review 
and further references). For a stellar cluster with density 
n* and velocity dispersion ofci — °i + °2 with two pop- 
ulations of stars with stellar masses mi and mi and stel- 
lar radii r\ and ri the co l lision time ca n be estimated as, 
ijBinnev fc Tremaind[l987l : lFreitagll2008t ) 



t co ii ~ 5 Gyr 



1 x 10 6 pc~ 3 



0~rcl 



2 R 2M P 







lOkmsec -1 ri + r2 mi + ni2 



(6) 



In Figure [7] we show this estimate for the collision time scale 
of putative star clusters forming when the gas would set- 
tle into rotational support in our simulations for our two 
assumptions regarding the loss of angular momentum. We 
take ZM/A-kR 3 as a proxy for the stellar mass density and 
estimate a va \ as V c /V2. For solar mass stars the collision 
timescales for a star cluster, with mass and radius similar to 
that expected for the gas in our simulations if it continues 
to loose angular momentum, could become shorter than the 
Hubble time (~ 4.16 x 10 s yr at z = 15). 

Stars in a populous cluster (iV 3> 10) exchange energy 
and angular momentum through two-body relaxation. The 
result of this relaxation process is the efficient mass segre- 
gation of the cluster with more massive stars concentrated 
at the centr e of the cluster and a much enhanced central 
density (e.g. ISpitzer|[l969l ; iBinnev fc Tremainej[l987l ) which 
eventually can lead to core collapse. The estimate in Figure 
[7] will thus be a substantial overestimate for the collision 
time of the central regions of such a cluster especially if core 

collapse has occurred. 

The relaxation time can be estimated as in iFrcitag 

iiooi), 



t rlx = 2Myr 



10 10 6 pc~ 3 
X nT 



10 km sec" 



1M 
< m > 



(7) 



If the gas in our haloes were to fragment and form a 



where A ~ Zn(0.02 N), N is the number of stars in the clus- 
ter and < m > is the mean stellar mass. The timescale for 
core collapse to occur is a few times the relaxation timescale 
at the centre of the cluster. If the timescale for core collapse 
is shorter than the main sequence life time of massive stars 
a very massive star should form by collisional runawa y col- 
lapse (|Giirkan et al.ll2006l ; IPortegies Zwart et al.ll2004h . 

In the bottom panel of Figure [7] we have plotted the 
relaxation time for our two assumptions of angular momen- 
tum loss. The central relaxation time scale and thus the 
time scale for core collapse of a putative cluster will again 
be substantially smaller than this. Conditions for the for- 
mation of a very massive star by collisional runaway and a 
subsequent collapse to a massive black hole appear thus to 
be favourable especially for s tellar clusters in the lower mass 
range around 10 4 - 10 5 M f 7 1 (|Freitadl2008l '). As discussed in 
£14.31 such a massive black hole would find itself in an ideal 
environment for rapid further growth. 
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Figure 8. Left-hand Panel: The fraction of enclosed mass as a function of the fraction of the virial radius. Right-hand Panel: The 
enclosed mass as a function of the fraction of the virial radius. The collapse into a compact configuration appears to be most efficient 
in haloes with V v i T ~ 20 km sec - 1 suggesting that massive seed black holes form most efficiently in these haloes. This should imprint a 
characteristic mass on the mass spectrum of massive seed black holes. 



4.5 Characteristic masses of massive black hole 
seeds 

So far we have concentrated on simulations A-C for DM 
haloes with virial velocities in the range 19 — 30 km sec" 1 . 
Recall that the mass of the discs in our simulations are 2 x 
10 4 Mq and 1 x 10 s Mq. Intriguingly the more massive discs 
form in the smaller mass haloes. In order to investigate if 
this is a general trend we have run two further simulations, 
D and E, for haloes with larger and smaller masses (see 
table[T]for details) extending the range of virial masses from 
~ 1 x 10 6 — 1 x 10 9 Mq which correspond to virial velocities ~ 
6—52 km sec -1 . Due to the adaptive nature of the Enzo code 
the spatial resolution and the gas resolution is comparable 
in all simulations. 

The fraction of the mass within a certain fraction of 
the virial radius at the time when rotational support sets is 
shown in the left hand panel of Figure [8] It appears to have 
a maximum for virial velocities around 20 km sec - . 

For lower mass haloes the temperatures attainable by 
atomic cooling is above the virial temperature of the halo 
and the gas cannot collapse without efficient metal or molec- 
ular hydrogen cooling. Surprisingly for the more massive 
haloes the collapse also becomes less efficient. 

In the right hand panel of Figure [8] we have plotted 
the total mass against the fraction of the virial radius. The 
collapse into a compact configuration appears indeed to be 
most efficient in DM haloes with virial velocities around 
20 km sec -1 and to lead to the formation of compact self- 
gravitating discs with characteristic masses of 10 M© . This 
should result in a peak in the mass function of black hole 
seeds and establish a characteristic mass for the population 
of massive black hole seeds. Predicting the space density 
of massive black hole seeds formed in the way discussed 
here would require to predict the fraction of haloes with 
V v ir ~ 20 km sec -1 for which the conditions for the suppres- 
sion of H2 and metal cooling are favourable which is beyond 



the scope of this paper. However a fraction as small as 0.1% 
would be sufficient to seed the population of present-day 
super-massive black holes in galactic bulges. 



5 DISCUSSION AND CONCLUSIONS 

We have discussed here the pathways to massive black holes 
in DM haloes with T v i r > 10000 K based on insights into the 
formation of self-gravitating compact discs from Enzo AMR 
simulations of atomic cooling driven collapse in such haloes. 

We argue that in a fraction of these haloes the condi- 
tions should be just right for a rapid and efficient build-up 
of massive black holes and rather different from those found 
in the DM host haloes of dwarf galaxies with similar virial 
velocities. 

For this to be the case the metallicity in the collaps- 
ing gas in these haloes has to be sufficiently low for metal 
cooling not to induce early efficient fragmentation and star 
formation. The gas can then lose 90% and more of its an- 
gular momentum due to supersonic turbulent motions be- 
fore it first settles into rotational support in a very compact 
configuration with rotational velocities well in excess of the 
virial velocity of the DM halo. The corresponding fat, self- 
gravitating discs are marginally stable against gravitational 
instabilities. 

If at this stage either H2 or metal cooling induce suffi- 
cient fragmentation/star formation a very compact nuclear 
star cluster would form offering favourable conditions for the 
formation of an intermediate mass black hole by collisional 
run-away possibly with the intermediate stage of a super- 
massive star. 

If instead after the onset of rotational support atomic 
cooling is still the dominant cooling process because the 
metallicity is still low and H2 cooling is suppressed due to 
external or more likely internal UV Lyman- Werner radiation 
the gas will still be unable to fragment while gravitational 
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instabilities will induce continued loss of angular momen- 
tum. 

If the entropy due to the released gravitational en- 
ergy can thereby be efficiently trapped the geometry of the 
inflow of gas will stay spherical and a supermassive star 
will form. The ongoing collapse of the outer mass shells 
occurs with mass accretion rates of a few tenth to a few 
Mp) yr~ . This i s just t he ran ge of accretion rate for which 
iBegelman et al.l (|2006l . 120081 ) have suggested that the su- 
permassive star should take the form of a rapidly accreting 
quasi-star with a pronounced core-envelope structure and a 
rapidly growing core which will turn into a rapidly growing 
stellar mass black hole. 

If the gravitational binding energy is instead trans- 
ported away in a disc-jet configuration the accretion flow 
driven by the in-fall of the outer mass shells is likely to be 
directed onto a remnant black hole of an ordinary star of 
which the DM halo should contain at least a few. 

In any case the infalling gas becomes strongly bound 
and has inflow rates which largely exceed the Eddington ac- 
cretion rate for a stellar mass black hole. Feedback due to 
individual supernova is expected to have little direct effect 
on the dynamical evolution of the gas at this stage and rapid 
growth to a massive black hole is expected. 

We have performed numerical simulations for a range of 
virial velocities from 6 to 52 km sec -1 . As expected in DM 
haloes with V V i T < 10 km sec -1 the temperature of the gas 
attained by atomic cooling only is too low for the gas to 
collapse. Somewhat surprisingly in the more massive haloes 
with masses above this threshold the settling of the gas at 
the centre of the DM haloes becomes rapidly less efficient 
with increasing mass. If the gas temperature is well be- 
low the virial temperature of the DM halo the gas appears 
to break more easily into several lumps and random mo- 
tions become more efficient in stabilizing against collapse. 
DM haloes with virial velocities close to the threshold of 
V c > 10 km sec -1 show the largest gas mass settling into a 
self-gravitating disc. 

This suggests that massive seed black holes will form 
with a characteristic mass. The exact value of this charac- 
teristic mass is uncertain but should be set either by the 
competition of accretion onto the massive black hole and 
consumption of gas in star formation or by the feedback on 
the inflow from the energy released in the accretion flow; 
10 5 — 10 6 Mq appears to be a likely mass range. The pop- 
ulation of massive seed black holes forming in pre-galactic 
haloes at z ~ 15 discussed here will naturally grow into the 
observed population in observed galactic bulges during the 
rapid hierarchical build-up of galaxies expected to occur at 
somewhat smaller redshift. The early formation of such mas- 
sive black holes strongly alleviates the many problems which 
have been identified if supermassive black holes had to grow 
from stellar mass black holes in smaller dark matter haloes 
or more generally in an environment where supernova feed- 
back from star formation can prevent the supply of sufficient 
amounts of fuel for rapid growth. 
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